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Abstract 

In this paper an active flutter suppression problem is studied for a thin airfoil in un- 
steady aerodynamics. The mathematical model of this system is infinite dimensional, 
because of Theodorsen’s function which is irrational. Several second order approxima- 
tions of Theodorsen’s function are compared. A finite dimensional model is obtained 
rom such an approximation. We use control techniques to find a robustly stabi- 
hzing controller for active flutter suppression. 


1 Introduction 

In this paper an active flutter suppression problem is studied for a thin airfoil in unsteady 

aerodynamics. Because of the interaction between the structure and the flow, flutter 
(dynamic instability) occurs at a certain flow speed. Therefore, it is important to de- 
sign active feedback controllers stabilizing the airfoil. A ribustly stabilizing feedback 
compensator is obtained from the H°° control theory. This theory gives us the largest 
amount of uncertainty (due to neglected aerodynamics) which can be tolerated in 
the problem of active flutter suppression. 

In general, mathematical models for airfoils in unsteady aerodynamics are linear 
time invariant infinite dimensiona l systems. The basic difficulty in such systems is 
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to compute the aerodynamic loads due to unsteady flow. The simplest models (in 
the frequency domain) for the unsteady aerodynamics contain Theodorsen s functl °" 
as the infinite dimensional part. There are several techniques for designing feedback 
controllers directly from the infinite dimensional airfoil model see e.g. [1J. In this 
method the controller itself is infinite dimensional, and hence one has to approximate 
it in order to obtain an implementable finite dimensional controller. Another method 
is to approximate the infinite dimensional part of the system and design a finite 
dimensional controller from the finite dimensional approximate model. In this paper 
we consider the second method, and design a robust controller, which stabilizes not 
only the finite dimensional model, but also the infinite dimensional model. The main 
tool used here in the robust controller design is the H°° control theory. 

In the next section we define a mathematical model for a thin airfoil. Several 
second order approximations for the Theodorsen’s function are compared in Section 
3 In Section 4 we present a robust stabilization algorithm for flutter suppression in 
the presence of unmodeled aerodynamics. Concluding remarks are made in the last 

section. 


2 A mathematical model for the airfoil 

We consider the following mathematical model (see e.g. [1], [2]), for a thin airfoil 
shown in Figure 1, 

M,z{t) + B.z{t) + K,z(t) = -J-fXO + ^ 


where z(t) = \h(t), a{t), 0{t)) T , and u{t) represents the control input. 
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where all the constants are related to the geometry and physical properties of the 
structure. 


In order to apply Laplace transform techniques, we will assume that z{t) = 0 for 
t < 0. This corresponds to the indicial problem (see e.g. [1], and [3]). Aeroelastic 
loads are represented by F(t) = [P(t), M a (t), M 0 (t)] T . We can represent F{t) as 


F{t) = M a z(t) + B a z(t) + I< a z(t) + F c {t) (2) 

where F c (t) is the “circulatory” part of F(t). The matrices M a , B a and I< a can be 
computed in terms of the problem data [9] [8] 
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where T,’s are Theodorsen’s constants, see e.g. [9], 

Using Theodorsen’s formulation, F c (t) can be expressed in the frequency domain 
as (see e.g. [9] pp. 395-396, or [8] pp. 26-28) 

F c (s) = C(s)(B c i + sB c2 )z(s ) /o\ 


C - d . 
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where 5 is the Laplace transform variable, ~ represents the Laplace transformed ver- 
sion of a time signal, C(ju) is the Theodorsen’s function, and are con- 

stant matrices given by B c i = b x c\ and B c2 = 6 ic 2 where 61 - pVb 


2nb(a + 0.5) 
T x2 b 


d = V [ 0 1 Tio/tt], and c 2 = [1 6(0.5 - a) bT u /2n}- 

Suppose that the measured output for feedback is 

y{t) := c-iz(t) 4- c 2 z(t). 

Then, taking the Laplace transforms of (1) and (2), and then using (3) we obtain a 
transfer function from u to y, denoted by P{s)- 

C 0 (sl - A)~ X B 


y( s ) 


= P{s) = 


G( s ) \-C 0 (sI -A)~'B l C(s) 

where C(s) is the Theodorsen’s function, and 
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Note that the plant can be seen as a feedback system whose feedback path consists 
of the aerodynamics represented by Theodorsen’s function, as shown below. 



Figure 2: Structure of the plant 

The function C(s) is irrational, and in practice it is approximated by a low order 
rational function, say C a (s). This leads to an approximate model for the plant to be 

controlled 

C,(sI-A)-'B, 

K{S) ~ C«(s)‘ 
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In Section 4 we will see that a rational feedback controller A'„, which stabilizes 
P a , stabilizes the original infinite dimensional model P if and only if the H°° norm 

of the “weighted” closed loop transfer function is less than the inverse of the L°° 
approximation error 


\\c - Calico := sup | C(juj) - c a (ju ) |. 

tjj 

Therefore, we have a better chance of stabilizing P by a rational K a if we can make 
the L ^ error in Theodorsen’s function approximation. In the next section we compare 
the L 00 errors of several second order approximates of the Theodorsen’s function. 


3 L Approximation of the Theodorsen’s func- 
tion 

As mentioned above, the Theodorsen’s function, C(s) which appears in the feedback 
path of the plant model, is infinite dimensional. For controller design (synthesis) and 
simulation (analysis) purposes we would like to use a finite dimensional approximate 
C a (s) instead of the exact irrational C(s), which is given by (see e.g. [9]) 

C{ju>) = Re[CO'w)] -I- jIm[C(ju>)] ^ 

where 

Re[C(;w)l = + jjM) + h(^)(Ti(^) ~ JqM) 

(Ti(o;) + J 0 (u )) 2 4- (Ti(u;) - T 0 (^)) 2 

Im[C(?'a;)l = (h(^)> oM + Ji{uj) JqM) 

(M”) + J 0 (w)) 2 + (Tj(w) - Jq(ui)) 2 

{ J o, J u Y 0 , Kj are the Bessel functions). Several second order approximations of (5) 

can be found in the literature, see for example [8], These approximations are in the 
form 


CJs) = + 

(1 + r 3 s)(l + r 4 s) 


( 6 ) 
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where ti ,t,, r 3 , r, are positive real constants to be chosen. For example the following 
sets of numerical values are proposed by R. Jones, W. P. Jones and R. L. Moore 
respectively 

Tl = 18.6 , r 2 = 1.97 , r 3 = 21.98, r 4 = 3.33 ( 7 ) 


n = 20.62 , t 2 = 1.85 , r 3 = 24.39 , r 4 - 3.125 
Tj = 10.61 , t 2 = 1.774 , r 3 = 13.51 , r 4 = 2.744 . 


( 8 ) 

( 9 ) 


For each of these sets of numbers the error function \C{joj) - C a {ju)\ is plotted in 
Figure 3. 



Figure 3: Error function for R. Jones, W. P. Jones and R. L. Moore 
approximates of the Theodorsen s function . 

As we can see from this figure, R. Jones ’s approximation is the best one (in the L°° 
norm) among the three second order approximates listed above. In different norms 
other approximations may be better than the one which is best in norm But 
since we are going to use H°° control techniques (in order to guarantee the robustness 
of the controllers derived from the approximate plant), we will need an error oun 
in the L°° norm. Below we will show that it is possible to improve the L°° error ot 
the R. Jones approximation by fine tuning the values of r x , . . . ,r 4 . 
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We found that the values of 


ri = 18.57 , r 2 = 2.057, r 3 = 21.93 , r 4 = 3.446 (10) 

give a function C a {ju>) whose magnitude is “close” to being a Chebyshev approxima- 
tion for the magnitude of C(ju>), (i.e. the error function | \C a {ju)\ - \C{ju>)\ |, shown 

in Figure 4, “nearly” satisfies the necessary and sufficient conditions for |C a 0w)| to 
be a Chebyshev approximation of \C(ju)\). 



Figure 4: 


\Ca(M\ - \C{jio)\ 


versus log(ai) 


We have obtained the above values for r.’s by slightly modifying the approximation 
scheme proposed in [11], We would like to determine if this choice for C a is a “good” 
L°° approximate of C. For this purpose we first point out the following relationship 
between the L 00 error and the error in magnitude and phase functions: 

Lemma: Let and 4> a (uj) denote the phase of C(ju) and C a {ju) respectively, i.e. 

C{ju) = \C{juj)\e^ C a (ju) = \C a {juj)\e^\ 

Then we have 


\C(ju>) - C a (ju>)\ < 


\C(M\ - \C a (ju)\ 


+ \C<jw)\ 


<j)(u) - 4> a {ui) 


( 11 ) 
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Proof: By definition following equalities hold 


| C(ju) - C 0 (ju,)| 


\C(ju)\e^ - \C a (ju)\e^ 

|C(ju>)|e i( * ( " )_ * o(u;)) - | C(ju)\ + \C{ju>)\ - |<?o(jw)| 

|C(;w)|(c jW( ‘ ,,) "*“ (w)) - 1) + (|CO w )l “ \C*U» 

On the other hand, for any 9 € [-* , A we have 
\e> 6 - 1| < |0|- 
Hence we conclude that 


| C(ju) - C a {ju >) I < 


\C(ju))\- \C a (ju)\ 


+ \C(ju)\ 


4>{u>) — 4>a{u) 


□ 


This lemma says that the Chebyshev approximation error 2 for the magnitude 
function plus the corresponding “normalized” phase error is an upper bound for the 
overall L°° error. We also deduce from this lemma that in order to ge a goo 
error bound we may try to develop an approximation scheme such that whenever the 
magnitude error is large, the normalized phase error is small and vice versa °wever, 
if we obtain C a from the Chebyshev approximate of \C{ju)l this automa ica \ 
the normalized phase error function, which does not necessarily satis y t e a o\ e 
mentioned nice property. However, we will see from the following numerical example 
that this property is satisfied for the second order approximation we have proposed by 
( 10V For C a determined from (10), the two terms in the right hand side of (11), as well 
as the function in the left hand side of (11), are shown in Figure 5. It is quite surprising 
that the normalized phase error function alternates with the magnitude error function, 

’What we mean by Chebyshev approximation for the magnitude function u i the ■ Sup- 

pose |C(jw)| is known, and we want to approximate the real valued function in the L norm b a 
function |C.(ju,)|; the problem is to find a real rational C.(.) (whose order is fixed) achieving 

smallest Chebyshev error 


sup 


|C(jw)| - |C«(j«)| 


In the text we use the term L« approximation for the approximation of the complex valued function 
C,tl 1 we - the term Chebyshev approx, motion for the approx, motion of the reo, valued 

function |C(jw)|. 
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i.e. whenever the first term is large the second term is small and vice versa. Also 
interesting is the fact that for this choice of C a the function \C(ju) - C a {ju)\ is an 
envelope of the two functions appearing in the right hand side of (11). 



the magnitude error and normalized phase error functions. 

The above observations can be generalized as follows: 

Conjecture: 

Let F(s) € H°° be a minimum phase and positive real function (possibly 
irrational), and let F a (s) <E H°° be a fixed order real rational function 
approximating F(s). If |F a (>u;)| is the best Chebyshev approximation of 
|T(ju;)|, then F a (s ) is the best L°° approximate of F(s). □ 


4 Active flutter suppression 

Let us consider the thin airfoil model obtained in Section 2. When flutter occurs the 
plant P(s) is unstable, and we would like to design a feedback controller stabilizing 
the closed loop system, shown in Figure 6. In our design we will use C a given by the 
numerical values in (10). This gives us an approximate plant model P a . A robustly 
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stabilizing finite dimensional controller K.(») will be obtained from P„, and it t will 
be shown that under a certain condition, this controller stabrhzes the ongmal infinite 
dimensional airfoil model, with a certain robustness level. 


P(s) 


Kjs) 


Figure 6". Feedback control system 


Consider the approximate plant 
C 0 {sl -A)' l B 0 


Pa(*) = 


1 — C 0 {sl “ A) 1 B\ C a (s) 

We can find rational transfer functions N U N 2 ,M G H°° such that 
C 0 (sl - A)~ l B 0 = ^7 and C °( sI ' A ^ Bl = ) 


M{s) 

Therefore we can express P and P a in the form 

N 0 (s) 


M{ s) 


Pa{s) = 


P(s) = 


N 0 (s) 


M{s) — N\(s)C a (s) ' M(s) — Ni{s)C{s) 

Thus, P and P a differ in their denominator, in the sense that 

NJs) , p , X _ Ap(s) 

P{S)= M^) a{S) " M pa (s) 

where N p {s) = N 0 {s ), M p {s) = M(s) - {s)C{s), M pa {s) = M P ( S ) + and 

& M (s) = N l (s)(C(s)-C a (s)). 
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Let t a be an upper bound of the L 00 approximation error for the Theodorsen’s func- 
tion, i.e. 


\\c - c a 


OO 


< Co, 


(note that for the choice of C a given by (10) we can choose e a = 0.012). Then from 
[4], [5] and [12] we C3,n deduce thcit & controller A ^ stabilizing Pa c^nd Achieving ctn 
H°° performance 


7 (A'„) = ||JV,M-'(1 + P.K.) 


-ll 


( 12 ) 


stabilizes the infinite dimensional plant P if 3 


1 


l(I<a) < -• 


(13) 


One proves this as follows: A controller K a stabilizes all plants of the form P a 


xr ~~r — if the roots of 

Mpa-L^M 


1 + K a (s)- 


N p (s) 


= 0 


M pa (s) - A m{s) 
are in the left half plane. This condition is satisfied if 


L a := 1 - + P a K a )~ l 


is invertible in H°°. Since HAjv/Hoo < e a |A r i( - /o;)| and K a stabilizes P a (meaning that 
(1 + P a K a )~ l € H cc '), a sufficient condition for L A to be invertible in H°° is (13). 

In fact, if K a stabilizes P a and satisfies (13), then it stabilizes P with a certain 
robustness level, see e.g. [7]. The controller I\° pt , which minimizes 7(A' a ) over all 
controllers stabilizing P a , has the best chance of satisfying (13). Note also that we 
increase our chances of satisfying (13) by decreasing e„. 

An interesting question about the stabilization of P by K a is: How much can we 
increase t a so that 


7 a •“ 


Ka 


inf „7(A'„)<- 

stabilizing r a 6 a 


? 


3 We would like to point out that the perturbation in the plant is in the denominator only, so the 
term y(K a ) is slightly different than the one in [5] and [12], where both numerator and denominator 
perturbations are considered. 
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The answer to this question gives the largest L°° error we can tolerate in approxi- 
mating C(ju>) so that we can still find an active feedback controller stabilizing the 
original plant. The problem of minimizing 7(A' a ) over controllers K a stabilizing P a 
is a special case of a one block H°° optimal control problem, and can be solved easily 
by finding the singular values and vectors of a Hankel whose symbol is a rational 
function, or by using the Nevanlinna Pick interpolation, see e.g. [4] and references 
therein. 


5 Concluding remarks 

An active controller design method is illustrated for a thin airfoil. The model P for 
this system is infinite dimensional. By approximating the infinite dimensional part 
of the plant we have obtained a finite dimensional approximate model P a . We have 
illustrated that using a Chebyshev approximation for the magnitude function \C(ju,)\ 
we obtain a finite dimensional approximate of C(joj) which is nearly optimal in the 
L°° norm. 

A finite dimensional controller K° pt can be obtained by solving the one block H 
problem posed in Section 4. In the H°° problem formulation we used the finite dimen- 
sional approximate model P a . We have shown that if the H°° optimal performance ; a 
is less than the inverse of the I 00 approximation error of the Theodorsen’s function, 
e a , then the controller I\ a stabilizes not only the finite dimensional model P a . but 
also the original infinite dimensional model P. 
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